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Abstract 

(N 

83m Kr is supposed to be used to study the properties of the windowless gaseous 
tritium source of the experiment KATRIN. In this work we deduce the amount 
of 83m Kr which is necessary to determine possible potential inhomogeneities via 
conversion-electron- line broadening. 
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vo : 1 Introduction 
o 

In the windowless gaseous tritium source (WGTS) for the experiment KATRIN, there 
should be - in the ideal homogeneous electrostatic potential which does not 

Ph! depend on time during the exposure in one spectrum point. In fact, this does not need 
Q_i! to be so. It might be both inhomogenous (due to, e.g., charging caused by the tritium 
decay or change of source's wall work function along the WGTS) and time dependent. In 
this work, however, we consider the inhomogeneity only. 



x. 

2 General 



To examine the inhomogeneity of the potential, we insert some source of monoenergetic 
electrons into the WGTS. In particular, we use gaseous 83m Kr. There exists a highly 
converted 2 isomeric E3 transition with the energy i? 7 =32151.5 eV The conversion 
electrons from the L3-subshell which are most suitable for our purpose have the kinetic 
energy i? =30473.1 eV (using the binding energy of the line width w=1.19 eV [5], 
and the corresponding relative conversion coefficient is 0^3/0^=0. 381 [I]. 

If the conversion electrons are being born in a region of inhomogeneous potential, 
their registered kinetic energies vary. This leads to line broadening and also to a line 
shift. Our task is to find the minimum quantity of 83m Kr needed to detect potential 
inhomogeneities at the level of ~50 mV. 
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2 a tot =2011 - evaluated by program [T] using the atomic potential of |2] 
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Let us imagine that the WGTS is cut into slices - as a salami - of the widths of dx. 
Then, at the position x, there is a potential distorsion 8jj{x) and the tritium gas density 
p(x). Let the probability that the electron undergoes i non-elastic scatterings prior to 
exit the WGTS be Pi(x). Taking into account the small width of the studied L 3 line and 
the finiteness of the electron energy loss in one scattering we can neglect all p^s except 
Po- 

The conversion line is described by a Lorentzian function. Then the count rate of 
the conversion electrons born in the interval (x, x+dx) is 

K(E;E ,w,x) dx = HK J ^\ J dx. 1 

It is easy to find that the correct normalizing constant, C, is given by 

c = A x (*L 3 /atot ^ 
Io L p{%) dx 

where A is the 83m Kr activity in the total volume of WGTS and integration runs over the 
whole length of WGTS. 

To achieve the total L 3 count rate, the formula must be integrated for x running 
over the whole WGTS length and we get 

V ; a to t 2tt Jo [E - E - e6u(x)} 2 + (w/2) 2 ' Jo HK 1 K J 

The KATRIN spectrometer is an integrating one. To express the shape which will be 
registered, the differential form must be convoluted with the instrumental response 
function j^j. This depends on the magnetic field intensities in the analyzing plane (Ba), at 
the source position (B s ), and the maximum intensity (B max ). Finally, we must multiply 
the result by the reduction factor due to the spectrometer acceptance angle |Sj 

Race = 7^(1 - cos$ max ) (4) 



where smtf max = JB S /B, 



3 Method 

In the experiment, we can detect the existence of inhomogeneities by measuring the kryp- 
ton line with and without, respectively, the presence of tritium in the WGTS and then 
mutually comparing the spectra. (No tritium implies no charging and then the inhomo- 
geneities are at least suppressed). Moreover, utilizing the fact that the detector is position 
sensitive (pixel-like), we can analogically search for different inhomogeneities in various 
parts of WGTS by mutually comparing spectra taken by different pixel groups. (Charac- 
terizing the inhomogeneities by some quantity AU - see below - we thus can distinguish 
if there are different AWs in different WGTS parts.) For the spectra comparison, we 
apply the statistical tests [7j. 
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To find in advance the sensitivity of this procedure, we apply the tests on simulated 
spectra corresponding to various kinds of inhomogeneities. Repeating such a procedure 
many times supplies a reasonable estimate of sensitivity. 

The shape of the disturbing potential, 5jj(x), is not known. Then it is necessary to 
use some reasonable estimate. We tried three different shapes (see Fig.flT])). 
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Figure 1: Shapes of disturbing potentials used in this work, (a) Triangular shape, (b) 
Gaussian one with e — 0.1, (c) Alternating Quenching Gaussian with e = 0.1. 



Triangular shape: the potential linearly decreases from the value in the tritium in- 
jection point (in the middle of WGTS), 5 v {L/2) = AU, to zero at both WGTS 
ends 

. , s / 2AU/L xx, < x < L/2 

du{x '-\2AU(l-x/L), L/2<x. (5) 



Gaussian shape: again, the potential decreases from AU but to some small value of 
^c/(0) = 8u(L) = e x AU at both ends and the shape is Gaussian instead of linear 
one 

5u(x) = AU x exp{-a(x - L/2) 2 }, (6) 
where a = -Alne/L 2 and its FWHM is 2^J\n2/a. 

Alternating Quenching Gaussian: this is Gaussian modulated by a periodical func- 
tion. It gives regions of small potential with opposite sign near the both WGTS 
ends 

S v (x) = AU x exp {-a(x - L/2) 2 } cos/3(x - L/2). (7) 

We request the minimum being (in absolute value) a small part of the maximum 
value, Su{x m i n ) = —e x AU and the potential to diminish at both ends. These 
two conditions lead to a system of transcendental equations for a, (3, and x m i n . 
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Fortunately, (3 = 3n/L may be extracted and a and x m i n may be separated. We get 
the equation 2 In cos sEmm _|_ f3x m i n x tan (5x m i n = which can be solved numerically 
and a = -f3/(2x min ) x taxi(3x min . 

4 Numerics 

The simulated spectra of the 83m Kr L3 conversion line were calculated in the energy 
interval from 30471 eV to 30475.4 eV with the step of 0.1 eV, i.e. in 45 points. The 
quantities characteristic for the experimental setup were taken those expected [S] to be 
used for KATRIN. In particular B max =6 T, B^=3 G, B s =3.6 T, and L=10 m. For the 
Gaussian and Alternating Quenching Gaussian shapes, e=10 _1 was chosen in both cases. 

The tritium density, p(x) is assumed to copy the tritium pressure dependence which 
we took following to Sharipov 9J. The probability Po(x) was calculated by Kaspar [TUj . 

To calculate the integrals in (jSJ), we utilized the Gauss- Legendre method as imple- 
mented in [TT]. The subsequent convolution with the instrumental function was done by 
simple summation using a suitably small step. The relative precision was checked to be 
better than 10 -3 which is quite sufficient for our purpose. The calculated spectra were 
then randomized in accordance with corresponding normal distribution. 

The tests were performed as follows: For a particular set of activity A, exposure 
time, and disturbing potential type eight spectra were generated corresponding to different 
A£7's. Then they were compared, each to each, using both \ 2 an d Student tests in both 
ratio and difference methods [7j. (Since all four tests gave comparable results we shall 
present only those of the \ 2 test of the ratio method.) The procedure was repeated ten 
times thus supplying a set of the test values - x 2 's per degree of freedom, xl, in this case. 
As the resulting value of "sensitivity" we adopt such difference AUi — AUj for which most 
of the xl values are greater than 1.5 and no one less than 1.3. 

5 Results 

We have considered all three types of the disturbing potentials described. In every case, 
the exposure time of 1 min in point, i.e. 45 min per spectrum, was assumed. The 
calculations were done for the 83m Kr activities of 10, 20, 50, 100, and 200 kBq, respectively, 
and counts registered by the whole detector were included. 

As mentioned above the spectra for every particular setup were simulated ten times, 
each comparison giving an interval 3 of the xl values. For an illustration, one such example 
is shown in Tabled In that case it is clear that the sensitivity of 40 mV is assured but, 
moreover, that of 35 mV can be expected safely enough. The sensitivities deduced in this 
way are presented in Fig. EJ Since we do not know exactly the shape of the disturbing 
potential it would be the most reasonable to use the worst (highest) values as estimates. 

It is necessary to remind here that the results are determined by the total number 
of counts registered. This means that the result for e.g. the activity A and exposure T 

3 E.g. for AUi - AUj = 0.02, we compare pairs of spectra for particular AU's (0,0.02), (0.01,0.03), 
(0.02,0.04) etc. We then present only the minimum and maximum of \t- 
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Table 1: Intervals of the xl values obtained during simulations. Activity od 20 kBq in 
WGTS, exposure time of 1 min in point, and detection by the whole detector are assumed. 



simulation 


AUi - AUj [V] 


number 


0.02 


0.03 


0.04 


1 


11-14 


11-19 


16-21 


2 


1.2 - 1.5 


1.3 - 1.9 


2.0 - 2.2 


3 


0.8 - 1.7 


1.7-1.9 


1.3 - 2.2 


4 


0.8 - 1.9 


1.1 - 2.1 


1.4 - 2.1 


5 


0.8 - 1.7 


1.0 - 2.0 


1.7-2.7 


6 


0.9 - 1.6 


1.3 - 2.0 


1.8 - 2.3 


7 


0.9 - 1.6 


1.3 - 2.1 


1.9 - 2.7 


8 


1.0 - 1.3 


1.1 - 1.9 


1.5 - 2.2 


9 


1.2 - 1.7 


1.4 - 2.1 


2.0 - 2.6 


10 


0.9 - 1.3 


1.0 - 1.8 


1.7-2.0 



are the same as those for the activity of A/2 and exposure of 2T. Similar argument holds 
when we compare the spectra measured by various parts of the detector. This is expressed 
in Fig. 121 by using at the horizontal axis the quantity 



», x, which is defined as 



x 



n . 4_ 
— x Ax -T, 

3 ' 



ntot 



where n to t is the total number of pixels in the detector, n is the number of active (regis- 
tering) pixels, A is the 83m Kr activity within WGTS in kBq, and T is measurement time 
(of the whole spectrum) in hours. 
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Figure 2: Sensitivities found for the three types of disturbing potential studied in this 
work. For the meaning of the variable x see text. 



5 



Acknoledgement 

This work was partly supported by the Grant Agency of the Czech Republic under contract 
No. 202/02/0889 and by the ASCR project AV0Z10480505 

References 

[1] M. Rysavy, O. Dragoun, M. Vinduska, Czech. J. Phys. B27 (1977) 538; 
O. Dragoun, M. Rysavy, J. Phys. G: Nucl. Part. Phys. 18 (1992) 1991. 

[2] C.C. Lu, T.A. Carlson, F.B. Malik, T.C. Tucker, C.W. Nestor,Jr., Atom. Data 3 
(1971) 1. 

[3] A. Picard, Z. Physik A342 (1992) 71. 

[4] K.D. Sevier, At. Data Nucl. Data Tables 24 (1979) 323. 

[5] J.L. Campbell, T. Papp, At. Data Nucl. Data Tables 77 (2002) 1. 

[6] A. Osipowicz et al. (KATRIN Collaboration), |arXiv:hep-ex/0109033l (2001). 

[7] O. Dragoun, M. Rysavy, A. Spalek, Nucl. Instr. Meth. A391 (1997) 345. 

[8] KATRIN Collaboration, "KATRIN Design Report", Report FZKA 7090, NPI Rez 
EXP-01/2005, MS-KP-0501. 

http : //www-ik . f zk . de/katr in/publications/ document s/FZKA7090 . pdf 
http : //bibliothek . f zk . de/zb/berichte/FZKA7090 . pdf 

[9] F. Sharipov, Numerical Calculation of Tritium Flow through the KATRIN Beam 
Line. 

http : //f isica.uf pr .br/sharipov/KATRIN/reportl .pdf 

[10] J. Kaspar, private communication. 

[11] M. Rysavy, Comp. Phys. Commun. 32 (1984) 191. 



6 



